Long wavelength helimagnetic order and skyrmion lattice phase in Cu20Se03 
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We report a long-wavelength helimagnetic superstructure in bulk samples of the ferrimagnetic 
insulator Cu2 0Se03. The magnetic phase diagram associated with the helimagnetic modulation 
inferred from small angle neutron scattering and magnetisation measurements includes a skyrmion 
lattice phase and is strongly reminiscent of MnSi, FeGe and Fei-^CoccSi, i.e., binary isostructural 
siblings of Cu2 0Se03 that order helimagnetically. The temperature dependence of the specific heat 
of Cu2 0Se03 is characteristic of nearly critical spin fluctuations at the helimagnetic transition. 
This provides putative evidence for effective spin currents as the origin of enhancements of the 
magneto-dielectric response instead of atomic displacements considered so far. 



PACS numbers: 75.25.-j, 75.85.+t, 75.50.Gg 

Major efforts have been made in recent years to un- 
ravel the nature of the magneto-electric coupling in mul- 
tiferroic materials [iHl]. Amongst a wide range of the- 
oretical scenarios two mechanisms are considered most 
prominent. First, a coupling mediated by effective spin 
currents in spin spiral magnets and, second, an exchange 
striction mechanism in which the magneto-elastic cou- 
pling proceeds via atomic displacements. An important 
property believed to provide unambiguous evidence of 
the latter mechanism is an enhancement of the magneto- 
dielectric response (MDR), describing changes in the di- 
electric polarization in applied magnetic fields or in the 
presence of long range magnetic order. However, the dis- 
covery of an enhanced MDR near the magnetic transi- 
tion of the insulator Cu2 0Se03 appears to question this 
view [5]. Detailed studies of the crystal structure and 
lattice dynamics strongly suggest the absence of sponta- 
neous lattices strains [SHH]. Being a lone pair containing 
piezoelectric ferrimagnet this was taken as evidence of 
a new magneto-electric coupling mechanism. Yet, en- 
hancements of the MDR without lattice strains are not 
specific to multiferroics and may represent a more gen- 
eral scientific challenge. For instance, the spin ice sys- 
tem Dy2Ti2 07 displays an enhanced MDR but is neither 
magnetically ordered nor multiferroic 0. 

Cu20Se03 is ideally suited to search for the origin of 
enhancements of the MDR without lattice strains. It 
crystallises in the non-centrosymmetric B20 structure, 
space group P2i3 [10], which structurally allows ferro- 
electricity. The unit cell is composed of three building 
blocks [llj. The first and the second building block are 
given by a square pyramidal and a trigonal bipyrami- 
dal CuOs unit in a three to one ratio, respectively. The 
third building block is a lone pair containing a tetrahedral 
SeOa. Magnetisation measurements and powder neu- 
tron diffraction have established ferrimagnetic order of 



the Cu^+ moments below Tc = 58.8 K, where three ferro- 
magnetically aligned Cu moments pair up antiferromag- 
netically with a fourth Cu moment O [12] . The exchange 
coupling is given by Jfm = — 50K and Jafm = 68K 
for the ferromagnetic and antiferromagnetic exchange, 
respectively, consistent with the Kanamori-Goodenough 
rules [T3] . 

Yet, the description in terms of ferrimagnetic order 
seems incomplete. First, well below Tc the magnetisa- 
tion increases almost linearly with increasing field be- 
fore saturating above ~ 120 mT without evidence for 
a spontaneous uniform magnetic moment at 5 = 
[3 [a US]. If the behaviour for B < 120 mT would 
be due to magnetic domains the slope and the onset of 
saturation would sensitively refiect demagnetising fields 
and depend on the sample shape which is not observed. 
Second, small changes of slope in the initial increase of 
M{B) suggest that the magnetic state is more complex 
[T4] . Third, the ferrimagnetic order is incompatible with 
the crystal structure of Cu20Se03 and a symmetry low- 
ering transition would be expected which is not observed 
[5] m [TS] m] . Finally, a detailed magnetic phase diagram 
of bulk samples was recently determined in magnetisa- 
tion and electric polarisation measurements for a (111) 
axis [T6[ [17] . Based on Lorentz force microscopy in thin 
Cu2 0Se03 samples helimagnetic order and a skyrmion 
lattice phase were attributed, but the phase diagram of 
the thin samples had completely different phase bound- 
aries. Moreover, for the thin samples the helical order 
was found to propagate along (110), which is not sup- 
ported by the B20 symmetry in bulk samples contrasting 
a key assumption in Ref . [TTl [18] . 

In this Letter we report a helimagnetic superstructure 
in bulk samples of Cu20Se03, which resolves all of the 
questions listed above and suggests strongly that the en- 
hanced MDR arises from spin currents due to nearly 
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critical helimagnetic spin fluctuations. The helimag- 
netic order in Cu20Se03 relates thereby to binary tran- 
sition metal (TM) compounds such as MnSi and FeGe, 
in which a hierarchy of three energy scales arises in the 
B20 crystal structure [T9|. On the strongest and sec- 
ond strongest scale these are ferromagnetic exchange and 
Dzyaloshinsky-Moriya interactions, respectively, gener- 
ating a long- wavelength helimagnetic modulation. The 
propagation direction of the helix is thereby the result 
of very weak magnetic anisotropics due to higher order 
spin-orbit coupling terms on the weakest scale. Most re- 
markable perhaps, a skyrmion lattice phase was recently 
discovered in the magnetic phase diagram of the binary 
TM B20 compounds ^2Qr.23j . which gives rise to a novel 
emergent electrodynamics [24l [23 . 

For our study a large Cu20Se03 single crystal (~ 
0.2 g) was grown by chemical vapour transport. Fur- 
ther properties of samples from the same growth were 
reported in Refs. |3 ^3 El ES] which also describe 
technical details of the growth procedure. An excel- 
lent sample quality was established by means of neu- 
tron diffraction at RESI at FRMII, notably an absence 
of impurity phases well below 1 % and a lattice mosaic 
spread smaller than the resolution limit of 0.02°. Using 
a neutron wavelength = 1.0408 A a simultaneous re- 
flnement of 60 Bragg spots provided lattice parameters 
a = b = c= (8.9199 ± 0.00104) A and a = /3 = 7 = 90°, 
in excellent agreement with the literature [TOl |TT] . 

The small angle neutron scattering (SANS) measure- 
ments were performed at the diffractometer MIRA2 at 
FRMII using a neutron wavelength = SA ± 5%. 
The neutron beam was collimated with an aperture, 
3 X 3mm^, placed '^1.5 m in front of the sample and a 
second aperture, 0.6 x 0.8 mm^, placed 0.38 m in front of 
the sample. The scattered neutrons were recorded with 
a CASCADE detector [27] located 2 m behind the sam- 
ple. For our measurements the sample was cooled with 
a Sumitomo pulse tube cooler, where a bespoke sam- 
ple stick permitted automated rotation with respect to 
an axis vertical to the incident neutron beam. The ro- 
tation axis coincided with a crystallographic (110) axis. 
The magnetisation was measured down to 3K with an 
Oxford Instruments vibrating sample magnetometer at 
magnetic flelds up to 9 T. The speciflc heat was measured 
in a Quantum Design PPMS using a quasi-adiabatic long- 
pulse technique. Laue x-ray diffraction was used to orient 
the sample. 

Shown in Figjl](a) are typical magnetisation data as a 
function of temperature in the vicinity of Tc. Well above 
Tc a strong Curie- Weiss dependence with a large fluctu- 
ating moment, /icw ^ 1.5/iB/Cu is observed in perfect 
agreement with the literature. With increasing applied 
magnetic fleld the magnetisation increases. In the vicin- 
ity of Tc faint maxima develop as illustrated in Fig{l](b), 
where M/B is shown for clarity. These features are anal- 
ogous to MnSi [28 , where they arise from the skyrmion 
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FIG. 1: (Colour online) Magnetisation of single crystal 
Cu2 0Se03 for various crystallographic directions, (a) Tem- 
perature dependence of the magnetisation in the vicinity of 
Tc. (b) Ratio /io M/B versus temperature revealing the fea- 
tures characteristic of the transition to the A-phase. Panels 
(c) through (j): Magnetisation as a function of field at vari- 
ous temperatures. Panels on the right hand side show typical 
data just below Tc, where a clear minimum in /io dM/dB, cal- 
culated from the magnetisation, is observed in the A-phase. 



lattice phase. The temperature dependence is consistent 
with the magnetic field dependence shown in Fig.[l](c) 
through (k) for field along (100), (110) and (111). With 
decreasing temperature M{B) increases before reaching a 
saturated moment = 0.48/iB/Cu at large fields. Ad- 
ditional changes of the slope of M{B) for temperatures 
just below Tc are best revealed in jnodM/dB^ where a 
distinct minimum is observed in a small T interval as 
illustrated in Fig.[l](e), (h) and (k). We thereby define 
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FIG. 2: (Colour online) Magnetic phase diagram of 
Cu20Se03 as a function of applied magnetic field for various 
orientations inferred from the magnetisation, (a) Overview 
for field parallel (111). Panels (b) through (d): Phase dia- 
gram in the vicinity of Tc for various orientations. Differences 
as a function of field are mostly due to demagnetising effects; 
the brown shading indicates the regime of nearly critical spin 
fluctuations. 



transition fields Bd^ Bai^ Ba2 and Bc2 as in the binary 
B20 compounds as shown in Fig.jl] (cf [29 ). 

Based on the magnetisation we obtain the magnetic 
phase diagram shown in Fig|2j The SANS data described 
below identify the following phases: (i) for B < Bd heli- 
magnetic order denoted 'h', (ii) for Bd < B < Bc2 con- 
ical order denoted 'c', (iii) for B > Bc2 field-polarized 
ferrimagnetic order, and finally (iv) a skyrmion lattice 
in the regime denoted 'A', just below Tc. We note that 
differences of Bc2 reflect demagnetising fields, which can- 
not be corrected accurately for the shape of our sam- 
ple. Likewise, the field range of the skyrmion lattice 
phase varies weakly with field direction (Figj2jb), (c) and 
(d)). However, the temperature range is clearly largest 
for (111) and smallest for (100), consistent with the mag- 
netic anisotropy favouring the propagation of the helical 
order at zero field along (100) [28l [3(H l3T] . 

Typical intensity patterns of integrated rocking scans 
that identify the various phases are shown in Fig. [3] Mag- 
netic rocking widths were small in all magnetic phases. 
An exception was the plane perpendicular to the applied 
field in the A-phase, where the precise intensity distri- 
bution was also sensitive to the field and temperature 
history. Future studies will have to establish whether 
this was the result of demagnetising fields related to 
the shape of our single crystal akin the shape depen- 
dence observed in MnSi single crystals [32 . For B = 
the intensity pattern consists of well defined spots at 
k - (0.0102 ± 0.0008) along all three (100) axes, 
characteristic of a modulation with a long wavelength 
Ah ~ 616 ± 45 A. This is shown in Figs.[3](a) and (b), 
which displays the intensity pattern for neutrons parallel 
(100) and (110), respectively. Preliminary tests with po- 
larised neutrons suggest a homochiral helical modulation. 
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FIG. 3: (Colour online) Typical integrated small angle neu- 
tron scattering rocking scans in Cu20Se03. (a) Zero field 
scattering pattern along (100) characteristic of helimagnetic 
order along (100). (b) Zero field scattering pattern along (110) 
characteristic of a helimagnetic modulation along (100). (c) 
Typical scattering pattern for Bd < B < Bc2 for T <^Tc. (d) 
Scattering pattern in the A-phase for magnetic field perpen- 
dicular to the neutron beam. Panels (e) through (h): Typical 
scattering pattern in the A-phase for magnetic field parallel 
to the neutron beam for various orientations. 



The very weak additional spots along the (110) axes in 
Fig.[3](a) are characteristic of double scattering. By anal- 
ogy with the binary B20 systems the scattering pattern 
at 5 = is characteristic of a multi-domain single- he- 
limagnetic state, where spots along each (100) axes cor- 
respond to different domain populations. In contrast, in 
MnSi the helical modulation is along (111). This is con- 
sistent with a change of sign of the leading order magnetic 
anisotropy in Cu20Se03 ["^ |2l] |30l |3T] , but contrasts 
distinctly the (110) propagation direction reported for 
thin samples [17]. 

In the range Bd < B < Bc2 the zero-field pattern 
(Figs.[3](a) and (b)) collapses into two spots parallel to 
the field. This is illustrated for B = 60 mT applied verti- 
cal to the neutron beam at 5 K in Fig.[3](c). Accordingly 
the modulation is parallel to B and, in analogy with the 
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FIG. 4: (Colour online) Typical temperature and magnetic 
field dependences of peak intensities recorded in SANS (see 
text for details), (a) Temperature dependence of the helical 
order, cf. Fig.[3](a). (b) Field dependence of the helical and 
conical state, where blue data points correspond the helical 
state, cf. Fig.[3](b), and green data points to the conical state, 
cf. Fig.[3](c). Labels (1), (2) and (3) mark the sequence m 
which sweeps were recorded; arrows show the sweep direction. 
Panels (c) & (d): Temperature and field dependence in the A- 
phase, cf. Fig[3](g). (e) Field dependence in the temperature 
range of the A-phase; helical order (blue, marked h), conical 
order (green, marked c) and A-phase (blue marked A), cf. 
Fig.[3](b), (c) and (d). 

binary B20 compounds, characteristic of a spin-flop phase 
also referred to as conical phase. 

In the A-phase, finally, the intensity pattern consists 
essentially of a ring of six spots perpendicular to the ap- 
plied magnetic field, regardless of the orientation of the 
single crystal with respect to the field direction. This is il- 
lustrated in Figs.[3](d) through (h). We begin with panel 
(d) which demonstrates, that the pattern for magnetic 
field perpendicular to the neutron beam is also perpen- 
dicular to the field. Further, Figs.[3](e) through (h) show 
the six-fold pattern for field parallel to the neutron beam. 
The orientation of the six-fold pattern in the plane per- 
pendicular to the field is thereby roughly aligned along 
(100), consistent with very weak magnetic anisotropy 
terms that are sixth order in spin-orbit coupling and the 
effects of demagnetising fields (see e.g. [20) [2T ). 

As demonstrated for the binary B20 compounds the 
six-fold pattern arises from a triple-/c state, with ki = 
0, coupled to the uniform magnetisation and stabilised 
by thermal Gaussian fluctuations (a single-^ modulation 
perpendicular to the applied field is energetically least 
favourable [20]). The topology of this peculiar multi-/c 
state is that of a skyrmion lattice, i.e. the winding num- 
ber is -1 per magnetic unit cell. This has been confirmed 
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FIG. 5: (Colour online) Specific heat, C, of CusOSeOs. With 
increasing temperature the helimagnetic transition is charac- 
terised by peak at Tc followed by a hump with a point of 
inflection at T2. The lack of field dependence at T2 is known 
as Vollhardt invar iance. The B and T dependence of C is 
characteristic of nearly critical helimagnetic spin fluctuations. 



experimentally in MnSi by means of Renninger scans in 
SANS [32 and the topological Hall signal pi. We there- 
fore interpret the A-phase in Cu20Se03 as a skyrmion 
lattice consistent with the study of thin samples [16j, 
where microscopic proof for the correct winding number 
is beyond the scope of our study. 

We also confirmed that the temperature and field range 
of the SANS patterns shown in Fig. [3] correspond with the 
phase diagram in Fig.|2j Typical T and B dependencies 
of peak intensities are shown in Fig. [4] (as peak intensi- 
ties are shown absolute intensities may not be compared 
easily in this plot). We find that: (i) the helical order at 
B = is characteristic of a second order phase transition 
at Tc (Fig.[4](a)); (ii) the transition from the helical to the 
conical state is at Bd and the suppression of the coni- 
cal state at Bc2 (cf. Fig.|4](b)); (iii) in the A-phase the 
signal of the conical phase vanishes completely in a nar- 
row range (Fig.[4](e)). It seems likely that the regime of 
the coexistence of the conical phase and the A-phase de- 
pends on demagnetising fields and thus the sample shape 
as recently observed in MnSi [28]). 

We return now to the origin of the enhanced MDR 
in Cu20Se03. At the accuracy of our SANS data the 
transition at Tc is second order (cf. Fig.[4|(a)). How- 
ever, the specific heat, C, shown in Fig.pFa) and (b), 
reveals that the transition consists of a narrow peak at 
Tc, characteristic of the latent heat of a first order tran- 
sition, and a broad hump with a point of inflection at 
T2. Under applied magnetic fields the peak and hump 
are suppressed with a shift of entropy towards high tem- 
peratures, while the point of inflection at T2 is invariant 
for B < Bc2 characteristic of a Vollhardt invariance at 
T2 [34 . For the case of MnSi it has been shown, that the 
same behaviour is due to a fluctuation-induced first order 
transition, where the chiral helimagnetic character of the 
fluctuations becomes dominant for T < T2. pQ) [35H38] . 
In turn, this suggests that the enhanced MDR arises from 
spin currents associated with the helimagnetic character 
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of the spin fluctuations. As the MDR is quantitatively 
rather smah, a fuh account connecting the magnetic with 
the dielectric susceptibility poses a challenge for future 
studies. Interestingly, the enhanced MDR in Dy2Ti207 
[9 quite likely originates in similar spin currents. How- 
ever, the chiral character of the spin fluctuations thereby 
must originate in geometric frustration rather than DM 
interactions. 

The remarkable consistency of the magnetic proper- 
ties of bulk samples of Cu20Se03 reported here with the 
binary B20 compounds is a surprise as the unit cell of 
Cu20Se03 contains 28 instead of 8 atoms. Thus, to the 
best of our knowledge, bulk samples of Cu20Se03 rep- 
resent the flrst example of helimagnetic order in a non- 
binary B20 compound, a B20 oxide, a B20 compound 
with a non-ferromagnetic leading order exchange interac- 
tion and a B20 insulator. The helimagnetism is thereby 
consistent with the proposal that thermal Gaussian fluc- 
tuations stabilise the skyrmion lattice In fact, our 
study clearly contrasts a theoretical proposal j39l |40| 
in which softened amplitude fluctuations stabilise the 
skyrmion lattice, since the softening would differ strongly 
between metals and insulators and with it the extent of 
the skyrmion lattice phase. Further, the helimagnetism 
resolves the open questions concerning the magnetisation 
of Cu20Se03, offering an explanation for enhancements 
of the MDR even for non-multiferroic systems such as 
Dy2Ti207. Perhaps most importantly, however, as an 
insulator the skyrmion lattice in Cu20Se03 promises an 
emergent electrodynamics akin that observed in its bi- 
nary metallic siblings [24, 25 , where electric flelds may 
now be used to manipulate the skyrmions [17 . 
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